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SECTION 1. INTRODUCTION 


The NMSU Center for Space Telemetering and Telecommunications 
systems is engaged in advanced communications systems research. Four 
areas of study that are being sponsered by this grant concern 
investigations into the use of trellis-coded modulation (TCM). In particular, 
two areas concentrate on carrier synchronization research in TCM M-ary 
phase shift key (MPSK) systems. A third research topic is the study of 
interference effects on TCM, while the fourth research area is in the field 
of concatenated TCM systems. 

Research on carrier synchronization in MPSK receiver structures has 
been motivated by the current interest in using 8 and 16 PSK TCM in 
power limited satellite channels to improve performance. NMSU has been 
involved in this research field since 1990 and prior to 1993 the main 
thrust of this research has concerned the theoretical analysis and 
simulation of MPSK carrier tracking. The carrier tracking approach that has 
been studied (and which is used in many BPSK and QPSK satellite 
receivers) is the high signal-to-noise ratio (SNR) approximation to the 
maximum a posteriori (MAP) solution to carrier phase estimation. The 
optimal MAP solution is a complicated mathematical calculation and was 
previously considered difficult to implement. With the advent of digital 
receiver designs this implementation barrier is now easily breached. This 
has motivated the NMSU Center for Space Telemetering and 
Telecommunications systems to fully analyze and study the optimal 
solution with the intent to simulate and implement it in hardware for 
study and evaluation at various SNR's persuant to the investigation of 
MPSK TCM systems. The second section of this report, "Optimal MAP 
estimation for MPSK carrier tracking" contains the results of this analysis 
which have been conducted to date. 

The hardware implementation of an MPSK carrier synchronization 
system has been in development at the NMSU Center for two years. The 
design has evolved out of a desire to conform to the above mentioned 

industry standard BPSK and QPSK carrier synchronization approach for 
comparison testing while being able to implement 8 and 16 PSK carrier 
synchronization. As a result the high SNR approximation was used for all 

four modulation formats in the first implementation of the MPSK carrier 

synchronization system. The system is under construction and should be 
completed by the end of July. Testing is scheduled immediately following 
construction completion and should be concluded in September. Carrier 
phase jitter and lock threshold measurements will be conducted in this 
first stage of testing. It is hoped that this data will verify previously 
collected theoretical and simulation data. The carrier synchronizer is 
described in section 3 of this report, entitled "MPSK carrier tracking 
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hardware project". It should be noted that the hardware implementation of 
the MPSK carrier synchronizer and the study of the optimal solution to 
MAP estimation are complementary. The digital design of the carrier 
synchronizer makes it a simple engineering change to implement the 
optimal solution and should time permit this change is scheduled for later 

this year. „ .. _ 

The fourth section of this report, "the effects of radio frequency 

interference on trellis coded modulation", concerns a simulation project. 
Radio frequency co-channel interference and burst interference are 
imparted upon a 16 state TCM system with Viterbi decoding and upon a 
Pragmatic TCM system. Using signal to effective interference ratios, bit 
error rates of the system were obtained as a performance measure, t is 
demonstrated that, in the presense of co-channel interference, 8 _ PSK TCM 
can achieve substantial coding gain as compared to uncoded QPSK. In the 
case of burst noise TCM systems perform badly. However with the help of 
periodic convolutional interleavers some of the losses can be recovered 

In the final section of the report, "concatenated trellis codes 
research work into combining TCM with Reed-Solomon (RS) codes is 
discussed. Previously, combinations of RS codes with binary convolutional 
codes have been used to achieve 10e-10 and 10e-20 BER s. The need for 
more spectrally efficiency has prompted investigations into exchanging the 
binary convolutional inner codes with Ungerboeck codes or the recently 
introduced pragmatic standard. This report describes simulation work and 
analysis being conducted to evaluate how well these concatenated systems 


perform. 
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SECTION 2. OPTIMAL MAP ESTIMATION FOR MPSK CARRIER 
TRACKING 

2.1 INTRODUCTION 

Technological advances continue to push the limits of the state-of-the-art in all areas. 
The area of communication systems is no different. As the demands of the “information 
age” require the transmission of greater amounts of data, the supply of resources is not 
necessarily able to keep up. This requires a more efficient use of the current facilities 
including satellites, transmission and reception stations, and existing equipment. 

Many current communication schemes use either binary phase-shift keying (BPSK) or 
quaternary phase-shift keying (QPSK) schemes for data transmission. While these have 
proven quite effective, the question naturally arises as to the possibility of increasing the 
data rate by applying higher-order phase-shift keying. Thus, an understanding of the 
higher-order phase-shift keying is essential in answering this question. 

The theory of Maximum A-Posteriori (MAP) estimation was developed in the ‘60s and 
‘70s.* At that time the implementation of the optimal estimators was not possible. As a 
result, simplifications were considered that facilitated the use of the basic principles of 
the MAP estimator, while admitting an implementation with available components. 
Thus, the full development of the theory was not pursued. Today, due to the advances 
in digital technology, the implementation of the optimal MAP estimator is now tractable. 
We thus consider how to implement the optimal estimator in a way that will readily 
retrofit onto current receiver structures. 

The chief result given here is a rederivation of the work of Simon using orthonormal 
basis functions that are associated with the mixers and integrate-and-dump circuits 
found in current, industry-standard receivers. The resulting error law that is derived is a 
function of the in-phase and quadrature signals that arise in a Costas crossover loop. 

2.2 MAP Estimation Derivation 

Consider the received MPSK signal, 


* The important works in this area are “Further Results on Optimum Receiver Structures for Digital Phase 
and Amplitude Modulated Signals,” (1978), M.K. Simon, Proceedings of the ICC, pp. 42.1.1-42.1.7, and 
“An Optimum Phase Reference Detector for Fully Modulated Phase Shift Keyed Signals,” (1969), S. Riter, 
IEEE Trans.. Aero, and Elect. Systems, Vol. AES-5, No. 4, pp. 627-631. 
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r(t) = cos (u> c t + 9 ( +9 m ) + w(t ) 


( 2 . 2 - 1 ) 


where E is the signal strength, u> c is the carrier frequency, 6 t is the phase shift in the 
communication channel, 6 m is the modulation angle, and w(t ) is additive white 
Gaussian noise with spectral density, N Q . It is assumed that the unknown phase shift, 
6j, is a random variable that is uniformly distributed between — 7 r and 7r radians. The 

receiver is thus assumed to have no prior knowledge of the channel phase shift. The 
receiver is to estimate 6 { by operating on the received signal r(t). The criterion of 

optimality chosen here, under which the estimate must operate, is the maximum a- 
posteriori (MAP) probability criterion. 


For the MPSK signal, in order to find the MAP estimator of the received channel phase, 
we must maximize p^|rj with respect to 9 i . By Bayes’ Rule, we obtain 


sup p 

e t 



max 

0 { € (— 7T, 7T ] 



max 

0 { e(- 7T,7r] 


p(r|*,)p(*,) 

p(r) 


( 2 . 2 - 2 ) 


Since the natural logarithm is a continuous, monotonically increasing function of a 
positive argument, then maximizing (2.2-2) with respect to 0 t is equivalent to 
maximizing the logarithm with respect to 9 r Thus, applying the logarithm gives 

1 M' , W r )]} = M' , ( r K)] + ln (' ,w < ) ) - ln ( p(r) )} (22 ' 3) 

The second term on the right hand side is constant, assuming that is uniformly 
distributed, and the last term is independent of 9 f . Thus, the maximization is really over 
the first term. The first-order necessary condition for optimality is 


<91n 



89 , 



= 0 . 


6<=e 0 


0i=0 o 


(2.2-4) 


In order to fulfill this necessary condition, we may work with the right hand expression, 
rather than the left. Hence, we must compute p(r|fy ). 

Recall that the received signal is given by 
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r (t)=lM. cos (u c t +e f +Q m ) + wit) 


(2.2-5) 


We may expand this signal using the orthonomal basis 

ifr 1 (t) = - ] pfcos(u c t+0 o ), ip 2 (t) = -^|r sin(u> c t +V- (2-2-6) 

r x = ^\ T Q r{r)cos{u c T + d 0 )dT, r 2 = ^ r(r)sin(u; c T + 0 o )dr (2.2-7) 

and obtain (in a mean-square sense) 

r(t) = r 1 ip 1 (t) + r 2 'i/j 2 (t). (2.2-8) 

Performing the integration for r x gives 

r i = ^fJ^-^^ C0S ( a; c^ + ^< J f 0 m ) + w{t) cos(w c f +0 o )dt 

= JjL J o r [cos(»„ + <« + C0s(2u«» + 9, + e o + e m )] dt + $ £.u*t ) cos(a; c ( + »„)<* 

= 41 cos(9„ + <P ) + -ft \l Ml) costal +e o )dl= s, + «j 
for u> c » 1 , where <f> = ^ — 0 O . A similar process yields for r 2 , 
r 2 = '\/^j ( ^['\/^ C0S ( CJ ^ + ^ + sin(w c f +6 0 )dt 

= - VFsin(0 m + <f>) + ^\*w(t)sin(u c t + 0 O )<* = 5 2 + n 2 . 

Thus, this expansion yields the “in-phase” and “quadrature” components of the received 
signal. These components are given with respect to the phase 0 O , which is generated as 

the output of a voltage-controlled oscillator (VCO). The in-phase and quadrature 
signals are currendy obtained in the industry-standard receivers. 

Since the noise was assumed to be Gaussian, and since ^{^i} = 5 i> and ^\ r 2i ~ 5 2’ and 
since the parameters r, and r 2 completely describe the received signal under the 

orthonormal expansion, then we obtain 
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(2.2-9) 


We have defined the variance of r(t) to be given by 

r 2 

a 2 = = +0 o )j > 

= rj T j & { w ( t ) u> ( T )} sin(u> c r-t-0 o )sin(u; c f + 0 o )drdt 

= - r) sin(u; c r + 0 O ) sin(u; c t + 0 O )drdt 

= ^ Jjsin 2 {u c t + 0 O )<* = 2S7 ( w c T " I sin(2w c T + 20 o )) 

= , w c »l- 

2 c 

Having found the probability density function in the general case, we may now consider 

each of the modulation shemes in order to evaluate the necessary condition. In each of 

the cases, we need to consider the probability density function of the modulation angle, 

9 in order to evaluate this necessary condition, 
m’ 

2.3 The BPSK likelihood Function Derivation 

For the case of binary phase-shift keying, we know that the modulation angle can take 
on 6 =0,7 r , and we obtain 

171 7 

p{e m ) = \8{9 m ) + \8{9 m -ir). 


Integrating with respect to the modulation gives 
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where 


X — (r 1 — $1 ) 2 + (r 2 - s 2 ) 2 = r\ + r\ + sf + s\ - 2r x s x - 2r 2 s 2 

=0 


0=0 


= r i +r 2 + E cos(p-r 2 sin<f>} = a-2^z 

where a = r\ + r\ + E , and z = r x coS(f) - r 2 sin^. Similarly, 
Y = {r l -s l ) 2 +{r 2 -s 2 ) 

= r\+r\ + E + 2^{r 1 cos<3 - r 2 sin < 3 ) = a + 2^¥z. 
Substituting X and yinto (2.3-1) yields 


fl=7T 


= f 2 +r 2 + s 2 + s 2 - 2r 1 s 1 - 2 r 2 s 2 


0— -TT 


■('M-Hii 


( -(<*-2 Jez)' 

(-(a+2-[Ez)'\ 


r^i r-v**y 

e- 2 *‘ ^ 

{ 2a 2 J 

+ e y ' 

-r 


l 


^4ttct j 



_ 


_ 


-&y 




2 o 2 J 


cosh 


-JI2 

2 

V <7 j 


(2.3-2) 


Defining now 

W = C0S( P ~ r 2 


= ^^jJ^^^r(t)cos(u; c / + 0 o }cos<f> - r(f)sin(u; c f +0 o )sin^)<# 
= At- ( t)c°s(u c t +0 o +(p)dt 

Y TtJ *0 

= j^^gr(t)cos(u c t +6 f )dt 


(2-3-3) 


where g = J^. We have introduced the constant £, since its presence in the form of 

an AGC circuit allows the resulting scale factors to be applicable to physical variables in 
industry-standard receivers. Substituting (2.3-3) into (2.3-2) gives 
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~(rl+r\ + E) 


± ±_ p "| 

p { r \ 9 i ,e o) = [j^Tj e 2 ° 2 cosh gr(t ) cos(u c t + O f )dt J = p{r\9 i j . (2.3-4) 

Not surprisingly, we see that ,0 o j is independent of 0 o . Thus, the first-order 

necessary condition becomes 


5? In ( r K) 


*i=*o 


d 

Mi 


-(f\ +r\ +£■) ^ 


In 

UO* 

f j \ 

+ lncosh| jpj Q gr(t)cos(u> c t + 9 { )dt j 



< J 



<>i=*o 


' T f* 

tanh gr(t)cos(u> c t +6 { )dt gr(t)sin{u c t +9^ dt 

. ™ o •'O y v o *0 


v i= v o 


2 E 
N„ 


tanh 


f 2 E ^ 
V N o , 


<3 = 0 


or IQ = 0 , where 


N. 


J = gf r(t)cos(<jJ c t +0 o )dt 

Q = -gj*r(t) sin (u c t + 0 o )dt = gjJVfr )cos(u c t +0 O +90 °)dt 


(2.3-5a) 

(2.3-5b) 


and /=tanh^^/J. A block diagram of this estimator may be obtained from this 
proceedure, and is shown in Figure 2.3-1 
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Figure 2.3-1 BPSK Carrier Phase Tracking Loop 


2.4 The QPSK Likelihood Function 


If the modulation is chosen to be quadrature phase-shift keying (QPSK) then we know 


that 6 m = or and 


p(# m ) = ?K - VJ + Mtn - 3 + Mm + VJ + S(6 m + 3 ^)]- (2-4-1) 


We now integrate with respect to the modulation probability and obtain 

f’( r K' 9 o) = l_,./'( r K' 9 o’ 9 ™) 

= = t) + = V-) + '>( r K = t) + 


m 4 




-X 


■fc)S' 


. 2a 


(2.4-2) 


where 


Xj = (rj - ) 2 + (r 2 - s 2 ) 2 


+ r 2 + Sj + s 2 ^- r \ s \ 2 r 2 s 2 

«n < 


<u=- 


" r i +r 2 + E ~r x 42E (cos^-sin<^) + r 2 V2£ r (cos^ + sind) 


(2.4-3) 


X 2 = (fj - s x ) + (r 2 - s 2 )" 


3 ^ - fj + r 2 + Sj + s 2 2/^Sj 2r 2 s 2 

{ m 4 


0 

m " 4 
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= r\ + r\ + E + r 1 V2£ r (cos^ + sin<p) + r 2 ^j2E (cosef) - sirup) (2 4-4) 

-^3 = ( r l — 5 l) +(r 2 —S 2 ) . _- 3 * =r \ + r 2 +S \ + S 2 ~ ^ r l 5 l ~ ^ r 2 y 2 _ 3 „ 

m 4 9 m=— 

= r i +r 2 + £ + r l 'j2E (cos(p-sin<p)-r 2 ^[2E(cos<p + sin<p) (2.4-5) 

*4 = (r, - ) 2 + (r 2 - 5 2 ) 2 = ^_ = rf + r 2 + sj + s 2 - 2r 1 5 1 - 2r 2 s 2 ^ 

m * ® m= T 

= r i +r 2 + £ (cos<£ + sin0)-r 2 V2F(cos0-sin<£) (2.4-6) 

Substituting (2. 4-3) -(2. 4-6) into (2.4-2) gives 

(2.4-7) 


(2.4-8a) 
(2.4-8b) 
(2.4-8c) 
(2.4-8d) 

(2.4-9) 

Under these definitions, we obtain 

W = 2cosh(r.t? 1 ) + 2cosh(oz' 2 ) = W x + W 2 . 

Simplifying, we obtain 


where we define 


a -r\ +r 2 + E , a = — and 

V 2(7 


2 Tj = -r^coscp - sincp) + r 2 (coscp + sin<p) 
z 2 = r l (cos(p + sin<3) + r 2 (cos<p - sin <p) 

2 t 3 = r^coscp - sin (p) - r 2 (cos<p + sin<3) = -z 1 
z 4 - -r^coscp + sincp)-r 2 (cos(p - sin <p) = -z 2 


az 1 az 2 az } az 4 
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W x - 2cosh|o(-r 1 (cos</> -sin^) + r 2 (sin^ + cos<^)jj 


= 2 cosh 


^^^|J o r r(0{ c °s(^ + 0 o )[cos <j> - sin<£] - sin (u c t + 6 0 )[cos4> + sin <£]}</? 


= 2coshj^-^#j J Q grit ){cos(u> £ / + 6 { ) - sin(w c / + 0, )}dt 
W 2 = 2cosh|a(r 1 (cos<£ + sin$) + r 2 (cos<fi - sin^))J 


(2.4-10) 


= 2 cosh 


r (t ){cos(o; c f + 0 o )[cos</> 4- sin^l)] + sin(u; c f +0 o )[cos<^ - sin <?!>]}<* 
J^r(0{ c °s(w^ + 0,.) + sin(u; c f + 0,)}<#j. (2.4-11) 


= 2 cosh 


Thus (2.4-7) can be written as: 

"(’M ) = (ifc 


f -a ^ 
V J 


(2cosh(£/) + 2cosh(\0) (2.4-12) 


where 


U = 


V = 



jj ^r(0{c°s(w c t + 0,)-sin(cJ c t + 6 i )}dt 
f grit )(cos(c*; c t + 6 ( ) + sin(u; t .t + 0,- )}dt 


We now consider the change of variables 


~u 


~X-Y~ 


1 -l" 

~X~ 

V 


X + Y 


1 1 

Y 


Solving for X and Y in terms of U and V , we have 


X~ 


1 -l‘ 

-1 

~U 

_/l\ 

‘ 1 l" 

~u 

- 

" u + v ‘ 

Y 


1 1 


V 

V 2 J 

-1 1 

V 

12J 

-u + v 


(2.4-13a) 


(2.4-13b) 


(2.4-14) 
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(2.4-15) 


gr (t ) cos(u> c f + 6 i )dt 
-jh Cs r ^ )sin(w c r + O i )dt 

V2<J 


We use the following property of the hyperbolic cosine 

cosh(jc — y) + cosh(a: + y) = 2cosh(.x)cosh(j') 

to see that 

cosh(f/) + cosh(V) = cosh(X - Y) + cosh(X + Y) = 2cosh(X)cosh(T). 
We now rewrite (2.4-12) as 




27T(7 2 J 


f ~(.r*+rl+E)^ 
e ^ 


(cosh(X)cosh(y)). 


Evaluating the first-order necessary condition gives 


dt l 



1 


-(r\+r\ + E) ^ 


H r M] 

II 

In 

U*> 2 ' 

+ In cosh(X) + In cosh(T ) 


l0 i=#o 


\ ) 





= tanh 


V* 

+ tanh 


— J-y j; gr (t ) cos (u c t + 9 ( )dt - -j3- J jjgr (t ) sin(u; c t +d i )dt 


J 0 i -Q o 


-^figriOsin^t+e^dt -ji-Y ^ gr(t)cos(u c t +G { )dt 

\ 


= ^ tanh 


( \ 


V^ 2 ) 


Q-^ytanh 


J6 i =6 i 


— E— o 


1 = 0. 


Thus, the necessary condition gives rise to the error signal 

tanh(^ lj<? - tanhf^Q j/ 


eyi) No 


where we may define 


(2.4-16) 


(2.4-17) 


-®o 


(2.4-18) 


PAGE 2-10 



Q = tanh 


I = tanh 


(2.4-19) 




A block diagram for the implementation of the necessary condition is given in Figure 
2.4-1. 



Figure 2.4-1. A Block Diagram Realization of the QPSK MAP estimator 


2.5 Higher-Order PSK 

We now consider the (higher order) M-ary phase-shift keying, where 

M = 2 n 2 < n - integer (2-5-1) 


and the modulation angle can take on the values 


d m =± 


( 2k + 1)7T 


k=0 y -,2 n 1 — 1, 


giving the probability density function 


2"- 1 -i 


k=0 v 


(2k+lh 




(2k+l)n 


Integrating with respect to the probability density function gives 

p[ r \ e i A ) = J_ p( r \ d i > 6 o> d m ) p( d m ) de m 
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2" *-l 


= iZ» &,AA 


k=Q 


(2&+1)tt 


.n-l ~*k 


= U _L_ y 

l 2" A 27t<t 2 ^ 


where for k = I,---, 2” \ 


2a y 

e + e 


(^) 


— (2fe+l)7T ^ 

2" J 


(2.5-4) 


k=\ 


X k — ( ^*i ^ r 2 $2 ) 


__ (24?+l)ir 


“ Jl 


2 2 2 2 ~ ~ 

= r i + r 2 + 5 i + 5 2 ” 2r i 5 i ~ 2r 2 5 2 


/i _(2fc+l)ir 
e ™~ 2 n 


= r 2 + r 2 + £ - 2r, VF cos(^ + *) + 2r 2 VF sin(2^ + * ) 
= r 2 +r 2 + S - 2r, VF j cos(^) - sin( (2fe ^ 1)7 ‘ ) sin(^) j 


+2r n 


, VF[ S in(a^)cos(^,) + cos(S^i) sin W) 
= a - 2VF(r ] [c t cos(^) - 5 fe sin(^)] - r 2 [5 (t cos(^) + C A sin(^)]j, (2.5-5) 


and where 


£ — 2 


n - 1 


2 2 

a = r, +r 7 + E , 


= cos(l^), S t 




(2.5-6) 


We also have for k — 1,- • •, £, 


— (r j Sj ) "t - (t~ 2 s 2 ) 


2 2 2 2 

= r, + r 2 + 5 l + $2 ~ 2r l 5 l ~ 2^ 2 5 2 


U _ (2*+l)» 
m ^ ^ 

= a- 2 4^[r 1 [c k cos (<f>) + S k sin(<£)] + r 2 [s* cos (<f>) - C k sin(<£)]) . 

We now define, for k = 1, •■•,£, 

z k = -rjc^ cos((p) - S k sin(^)J + r 2 \s k cos(^) + C k sin(^>)] 

z_ k = -rjc* cos (4>) + S k sin(0)] - r 2 [5* cos (<j>) - C k sin(<£)] 


(2 k + l)ir 


(2.5-7) 


(2.5-8a) 

(2.5-8b) 
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and note that for k = !,•••, i. 


X k = a- 2 -(e z k , X_ k = a- 2 41 z_ k 

Since Ci—k+i = ^t—k+i ^k j 

z k = ~ z -U-k+\Y 

Define now, 

W = £<■**+ A* 

k=\ 

and note that 

= Y j e PZk +/ JL(/ -* +1) = +e Pz - (e - k+1) + e fk ~ k +e* l ~ 


(2-5-9) 


( 2 . 5 - 10 ) 


W 


k=l 


k+1 


l A 


(2.5-11) 


= ^2cosh(/2r^) + 2cosh(j0z_^). 

k=i 

Examining the term 

pz k = -{3r^C k cos(<p ) - S k sin(<3)] + (3r 2 [s k cos(<£) + C k sin(^)] 

= -^-^^J o r r(r)cos(u; c T + 0 o )dr cos (</>) - S k sin(^)j 

+ r(r)sin(cj c r + 0 o )dr^S k cos(^) + C^sinG?)] 

= f T r(r){-cos(a; c r + 0 Q )\C k cos(<f>) - S k sin(<3)] + 
a •'O l L 

+sin(o) c r + 0 o )j^5^ cos (<p) + C k sin(<j!>) jj^/r 
= ^ [ T r(r)\-C k cos (u c t + 0 t ) + S k sin(u; c r + O^dr 

= -(¥)/'-(¥)<?'= - 7 *r - hQ' (2 5 ' 12) 

for 
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(2.5-13a) 


C b E 2 C b E S k l 

7* = -V = -TT-. H = — 

a Jy 0 o 


_S k E 2 S k E 


1 N 


I' = g]\(t)cos(u> c t +6 i )dt , Q' = )sin(w c f + 9 { )dt. (2.5-13b) 

Similarly, 

(3z_ k =-(3r 1 [c k cos (<f>) + S k sin(0)] - (3r 2 [s k cos(<j > ) - C k sin(0)] 


c,£ )/' + (¥)<?' = -V' + 


(2.5-14) 


Thus, W becomes 


l A A 

W = ]£2cosh(7^/' + n k Q') + 2cosh(-7 k I' + n k Q') = 4cosh(7*/0cosh(p A <2')- (2.5-15) 

k=i k=\ 


In this case, we may write the probability density function as 

„«-i f-*V| 


+c ^ i" 7 (2 - 5 - l6a) 


ln[p(r|^,0 o )] = ln (p^r } 
Evaluating the first-order necessary condition, gives 


Inh^-rV’' + H W Y 

V 2 7 TO J 


d In 

dW/ 

_ M 

Mi 

W 


0 l=°o 


(2. 5-1 6b) 


= 4 £ 7* sinh( 7/fe r ) cosh(/i*<2')f Jf] + cosh(7 Jfe / / ) sinh(/x*Q')f 

i i V y v 


In order to evaluate this expression, we note that 
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= -/ 


(2.5-17) 


f 


V 


dr 


A 


89. 


« J 





f 8Q'^ 

~d¥ i 


Q{-®o 



Oi~®o 
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Thus the necessary condition of optimaltiy gives rise to the error function 

K \ 

£ cosh( 7 */) cosh(;x ib Q)( 7jb <3 tanh ( 7 */) - fi k I tanh (/x*<?)) 

£<1) = - % 

£ cosh( 7fe /) coshf/x^Q) 


k=i 


or 


& tanh( 7 | ,I) - /q./tanM/q.Q) 

y cosh(jjI)cosh(fj.jQ) 

1 + ^ 1 cosh( 7jfe /)cosh(/x Jfe <3) 
j*k 


(2.5-19) 


2.6 A Comparison of the Various Strategies 

In comparing the various strategies, we observe the error function. In each case, it is a 
function of the in-phase and quadrature components of the signal. For higher order 
strategies, this becomes increasingly complex. We have primarily considered only 
strategies that are powers of 2, although others could be considered. This restriction 
corresponds to the sending of n bits of data in each symbol. This also gives rise to a 
great deal of symmetry in the derivation. 
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In Figures 2.6-1 through 2.6-4, we see the error signals for BPSK, QPSK, 8PSK and 
16PSK. assuming a 5 dB signal-to-noise ratio. In each case we see rounded edges 
associated with a “soft-decision” structure. Indeed, as the order increases, the edges 
seem to become more rounded, indicating the difficulty of making a decision. 

Another feature of the increasing order is that the magnitude of the error function 
decreases. Even though / and Q are plotted on the same scale, we see that the error has 
decreased by an entire order of magnitude from BPSK to 16PSK. This is also indicative 
of the conservatism required in executing the higher-order estimation. 
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Figure 2.6-4. 16PSK Error Function 
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SECTION 3. MPSK CARRIER TRACKING HARDWARE PROJECT 


3.1 INTRODUCTION 

To facilitate the testing of carrier synchronization techniques 
for MPSK systems a testbed is being constructed. The testbed is made 
up of two parts. The first of these is the carrier synchronizer. The 
second is the demodulator. 

The carrier synchronizer is capable of tracking BPSK, QPSK, 
8PSK, and 16PSK signals. The IF frequency is 15 MHz and the data 
rate can be 500k symbol-per second (sps), lMsps, 2Msps, or 4 Msps. 
The current synchronizer configuration uses the high SNR 
approximation to MAP estimation of carrier phase to track the 
received signal. It is a digital design. Its inputs are digital samples of 
the I and Q channels which it gets from the demodulator. This is 
shown in Figure 3.1-1. The synchronizer generates a signal which is 
transferred to a Stanford Telecom’s 1272 NCO board. This circuit 
generates the sine and cosine signals which are used to drive the 
demodulator. 

In Figure 3.1-2 the demodulator is shown. The input signal for 
testing carrier tracking is simply a 15 MHz carrier. Since an M-ARY 
PSK signal has M ambiguous lock points in its constellation 
transmitting the same symbol repeatedly is equivalent to randomly 
switching between all of the M symbols, in so far as carrier tracking 
is concerns. This makes it possible to create a simplified test 
environment utilizing nothing more than a 15 MHz carrier as an 
input. The demodulator utilizes the quadrature structure and both 
the I and Q branches are sampled to provide the digital information 
required for the carrier synchronizer. Both the demodulator and the 
carrier synchronizer require a symbol clock and this must be 
provided by a function generator. The system has been designed to 
allow for easy conversion to an actual MPSK receiver environment. 

3.2 THE CARRIER SYNCHRONIZER 

There are four digital circuit cards arranged along a backplane 
that make up the main part of the carrier synchronizer. These are 
shown in Figure 3.2-1. Referring also to Figures 3.1-1 and 3.1-2 there 
are several different functions provided in this unit. Starting at the 8 
bit I and Q inputs, two 512k EPROMs are used to generate the sine of 
the phase error during a symbol period. The 16 bit sine of the phase 
error is used as the tracking error signal and is passed to the loop 
filter. The loop filter’s transfer function is 1 + a/s and is implemented 
with the use of two 16 bit ALU’s; one for the accumulator and one for 
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the sumer. DC loop gain and the filter cut-off frequency are adjusted 
by the use of barrel-shifters in the loop filter circuit which can 
divide a digital word by powers of two. The input to the NCO is a 32 
bit word representing the frequency it is to output during the next 
symbol period. The 32 bit word must contain the quiescent 
frequency and any error signal component that the loop generates. 
The quiescent frequency is generated by the sweep circuitry and 
added to the error signal using two cascaded 16 ALU’s. The sweep 
circuitry can output a static 32 word for a quiescent 15 MHz or it can 
be swept to better match the test input signal. The sweep speed and 
step size can be changed from the front panel. Sweep controls are 
manual in the carrier synchronizer's current design. 

Figure 3.1-1 shows a phase shifter incorporated with the sweep 
circuitry. This was done to facilitate ambiguity resolution. The phase 
shifter can introduce a phase shift to realign the NCO frequency with 
a different MPSK phase vector. Since ambiguity resolution is not 
required for carrier synchronization the current synchronizer 
configuration does not include the phase shifter circuitry. However, 
its design has been completed for future use. 

The 32 bit word that is sent to the Stanford Telecom NCO board 
is transmitted in byte-serial fashion over an 8 bit bus. The four bytes 
are clocked into the NCO during a symbol period by the use of special 
timing circuits. These circuits use tapped delay chips and XOR gates 
to generate four clock pulses from a single symbol clock pulse. 

One other circuit that is incorporated in the synchronizer 
design is the lock detector circuit. In a future version of the 
synchronizer, lock detection will be automatic and tied into the 
sweep circuitry to initiate and terminate carrier sweeping. In a full- 
up receiver design this will assist in achieving symbol 
synchronization and ultimately in achieving carrier synchronization. 

The carrier synchronizer is a 10 inch high 19 inch wide rack 
mountable unit. The front panel design is shown in Figure 3.2-2. All 
of the user interface hardware (displays, switches, etc.), power 
supplies, fans, the card cage, and the 1272 NCO card have been 
installed in the unit. A 50 MHz oscillator source for the 1272 NCO 
card has also been included. With the exception of the phase shifter 
and lock detector circuits, all other circuits are completed. 
Preliminary open loop testing has been conducted on the 
synchronizer. This was done by replacing the EPROMs that contain 
the sine of the error function with EPROMs that contain a 
programmed sine wave. The sine wave is clocked into the loop filter 
and then to the NCO by means of a test counter circuit that’ is 
connected to the EPROM address lines. Using a 15 MHz quiescent 
frequency the carrier from the synchronizer swept back and forth 
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following the sine wave input to the loop filter. The next stage of 
testing will incorporate the demodulator and close the loop. 

3.3 THE DEMODULATOR 

The demodulator has both analog and digital circuitry. The 
front end is analog and comprised of off-the shelf RF circuit 
components. They have been selected and specified to meet power 
requirements in accordance with the input test signal, the sine and 
cosine from the carrier synchronizer, and with the inputs to the 
demodulator analog-to-digital converters . The test "received" signal 
is a composite signal constructed from the outputs of a 15 MHz 
frequency synthesizer and a NoiseCom white noise source (100 Hz to 
500 MHz). The signal and noise are summed and amplified. This 
composite signal is then split, 3 ways, to supply the two quadrature 
mixers and an output test port from the demodulator. The NCO sine 
and cosine outputs are separately amplified using two-stage fixed 
and adjustable amplifiers. They are then mixed with the composite 
"received" signal. The mixer IF outputs are low-pass filtered at a 2 
MHz bandwidth before entering the ADC circuits. The ADCs are multi- 
stage circuits and are the most complicated section of the 
demodulator. Both gain and DC offset must be achieved to drive the 
converters with the analog I and Q channel samples while 
maintaining amplifier stability. There are three op-amps in each ADC 
circuit (one for I and one for Q) to achieve the proper signal 
conditioning. The ADC's themselves, AD4098's, are clocked at 500 
kHz, the same "symbol" clock that drives the carrier synchronizer. 
The 8-bit ADC outputs are sent to the carrier synchronizer and 
represent the I and Q samples that are used as addresses to the 
EPROM's. The construction of the demodulator will be completed by 
the end of July and. 

3.4 TESTING 

Once construction, checkout, and calibration of both the carrier 
synchronizer and the demodulator are completed, MPSK carrier 
tracking tests will begin. The theoretical work and simulation results 
on MPSK carrier tracking, completed under NAG5-1392 will be 
verified by conducting phase jitter measurements. The signal-to- 
noise ratio (SNR) can be varied by attenuating the noise source in the 
demodulator. An HP 5172A time and frequency analyzer will be 
used to conduct jitter measurements. Further the lock thresholds of 8 
PSK and 16 PSK can be compared with those of the simulation 
results. 
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3.5 CONCLUSIONS 


After the phase jitter and lock threshold testing has been 
completed, it will be possible to further explore the use of optimal 
MAP estimation in MPSK carrier tracking as discussed in Section 2 of 
this report. By replacing the synchronizer EPROM's, which contain the 
high SNR approximation to MAP estimation, with EPROM's that 
contain the optimal solution, it will be possible to study the effects on 
phase jitter and lock thresholds that using the current standard of 
the high SNR approximation has on 8 PSK and 16 PSK. 

The entire system has been constructed with the idea that it 
can be used as an MPSK receiver with minimal engineering changes. 
Incorporating a phase shifter and lock detector with an external 
symbol synchronizer will complete the design of a fully functional 
MPSK receiver. This will be a valuable laboratory tool in the analysis 
of communications systems that utilize MPSK. 
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• Figure 3.1-1. The carrier tracking test system, emphasizing the carrier synchronization circuitry 
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Figure 3.1-2. The carrier tracking test system, emphasizing the demodulation circuitry. 





Figure 3.2-1 . MPSK Carrier Synchronizer Block Diagram. 
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Figure 3.2-2. MPSK carrier synchronizer front panel artwork. 










SECTION 4: THE EFFECTS OF RADIO FREQUENCY INTERFERENCE ON 
TRELLIS CODED MODULATION 


4.1 OVERVIEW 

Trellis coded modulation (TCM), pioneered by Ungerboeck 1* 2 has 
been cited as the modulation technique for the next generation of space 
communication systems. This being the case, the performance of TCM 
systems exposed to additive white Gaussian noise (AWGN) has been 
researched thoroughly. However, questions have surfaced concerning how 
well TCM systems will respond to radio frequency interference (RFI). 

This paper presents results of simulations geared to measure the 
performance of TCM systems in radio-frequency (RF) co-channel 
interference and burst interference. The software package. Block Oriented 
Simulation Software (BOSS), was used to simulate (at rate 2/3) a 16 State 
TCM system with Viterbi decoding and a Pragmatic TCM system 4. The 
performance of the systems is measured by their bit error rate (BER) as a 
function of signal to effective interference ratio. 

Results show that when compared to the uncoded, quadrature phase 
shift keying (QPSK) system, TCM achieves a substantial coding gain when 
subjected to a mixture of AWGN and co-channel interference. It is a fact 
that BER performance of TCM systems degrades significantly as a result of 
burst interference. However, periodic convolutional interleavers (PCI) can 
recover some of the losses. 

4.2 INTRODUCTION 

It has been effectively argued that Trellis coded modulation (TCM) 
provides a variety of advantages for space communications 5 & 6 j n suc h 
an environment (bandlimited, power limited and possible low signal to 
noise ratio (SNR)), TCM provides coding gain, and error correction without 
an increase in bandwidth. The bit error rate performance of TCM systems 
in Additive White Gaussian Noise (AWGN) has been well documented 
theoretically and by simulations 1, 2, 3, 5, 6, & 9. While trellis codes have 
been proven for Gaussian channels, their performance in non-Gaussian 
noise has been brought into question. In this project two cases of non 
Gaussian interference will be studied, co-channel interference and burst 
interference. 
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Co-channel interference is a typical source of degradation and has 
been the subject of some other investigations. For example, Losquadro et al 
published results on a study of RFI in DRS links and the impact on coding 
schemes used in that system 13. Other theoretical works have examined 
M-ary PSK subjected to AWGN and RFI 7. 

The second case to be studied is that in which the interfering signal 
occurs in bursts. It is already known that TCM performs poorly when the 
interference occurs in bursts ®. A popular scheme used to combat burst 
errors is interleaving. In a paper published in 1971, G. D. Forney detailed 
the structure of the convolutional interleaver used to combat bursts on the 
channel 10. A report published by Richer demonstrated improved 
performance with block interleaving for a Viterbi decoder subjected to 
burst interference 1 1 . Evaluations of how well communication systems 
work when subjected to various kinds of interference is an ongoing task. 
As new systems are developed, like TCM, they are put through the 'rings' 
to determine their limitations. 

The Telemetry Center at New Mexico State University (NMSU) has 
developed a data-base of TCM systems and simulation results in BOSS. 
From the data base two systems, 16-State TCM and Pragmatic TCM, were 
selected for this study. The details of their design can be found in ^ & 6 
Both systems make use of an 8-PSK (phase shift keying) signal set and 
their co-ordinates are transmitted over T and 'Q' channels ( in-phase and 
quadrature). Use of such a signal set for simulations lessens the amount of 
memory used, and the sampling requirements for signals. The TCM 
systems were re-configured to include a RFI vector source (co-channel 
interference model) and their bit error rate (BER) performance tested. 

Regardless of how the coded system is affected by various kinds of 
interference, it should be kept in mind that the value of coding is 
determined by comparing the performance of an uncoded system 
subjected to identical conditions. 

4.3. CO-CHANNEL INTERFERENCE 

To a user, co-channel interference is any unwanted signal energy 
within a frequency band of interest. The sources of such interference are 
many and varied but may be typified by signals from adjacent frequency 
bands. It is clear that whatever the source, such signals will cause 
degradation in the performance of communication systems. 
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In order to characterize the performance of TCM systems, the 
approach taken in this effort is to generate a parametric family of BER 
curves as a function signal to composite interference ratio. The total 
interfering power is the sum of AWGN and the RFI source in the channel. 
Each curve is classified based on the percent RFI energy content of the 
total interference. Comparisons are made with the corresponding uncoded 
QPSK system under the same conditions. For this purpose the systems were 
constructed in BOSS for simulations. In these BOSS systems signals are 
represented as vectors and thus create a demand for a RFI vector space 
model (see Figure 4.3-1). 


RFI Source 



Figure 4.3-1 TCM System with RFI Model 


4.4 RADIO FREQUENCY INTERFERENCE VECTOR MODEL 

Figure 4.4-1 illustrates the vector space of a QPSK signal seL A vector 
space model of RFI may be described as follows for co-channel 
interference. The interfering signal is considered as a sinusoid with fixed 
amplitude and fixed frequency within the transmission bandwidth, 

A*sin[ Wit + 0(t) ]. The phase is a random variable uniform from -k to k. 
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QPSK Signal Set 


Figure 4.4-1 RFI Vector Model 


The received signal is therefore the sum of the transmitted symbol, AWGN, 


and RFI. Correlation with reference signals '\j j. cos[ 2k f s t ] and 


Vf" s 


sin[ 2k f s t ], produces in-phase 


(I) and quadrature (Q) components of 


each term in the received signal. The correlator’s response to the RFI 
produces an I-term 


A sin[ 2k (fj - f s ) T + 6(t) ]- sin[ 6 ] 
T 2k (f| - f s ) 

and a Q-term 

A cos[ 2k (fj - f s ) T + 9 ] - sin[ 6 ] 
T 2k (fj - f s ) 


Adding the sum of the square of the above terms gives the energy which is 
reduced to A 2 sinc 2 [(fj - f s )T]. From this derivation it is justifiable to 
consider the effective amplitude of the interfering signal, the square root 
of the energy, as the magnitude of the above vector. Our model of the RFI 
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vector is therefore Aeff ei 0 . Consequently the model results in a complex 
vector. The magnitude of this vector is a function only of the difference in 
frequency (between the interfering signal and the reference), e is a 
random variable uniform from -u to n. With this complex vector T and 'Q! 
components are generated by extracting the real and imaginary parts. 

4.5 GENERATING THE RFI MODEL IN BOSS 

A block diagram of the BOSS modules used are show in Figure 4.5-1. 
The phase is derived from a random number generator whose output is 
uniform from -it to it. The phase term is made complex and is then 
multiplied by, ’Aeff', the magnitude of the interfering vector. 


Input 



QrComponent 


Figure 4.5-1 RFI Source Model 


Outputs from the 'Real' and 'Imaginary' modules, I and Q, respectively, are 
the additive interference to I and Q_ components of the message. With this 
new channel module inserted in the already built systems, their BER 
performances were measured. Plots of BER verses signal energy to noise 
ratio are shown in Figures 4.6-1, 4.6-2 and 4.6-3. The noise source is the 
sum of AWGN and RFI. 

4.6 CO-CHANNEL RESULTS 

Figures 4.6-1 to 4.6-3 show the parametric BER plots for the uncoded 
QPSK system and Trellis Coded 8-PSK systems. The TCM systems used are 
the 16-State system introduced by Ungerboeck 1 & 2 , and the Pragmatic 
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TCM, more recently introduced by Viterbi 4 . Plots are shown for RFI 
comprising from zero to one hundred percent of the total interfering signal 
in increments of twenty five percent. It can be seen from the plots that 
there are significant coding gains relative to the uncoded case. As an 
example, at 25% RFI for a BER of 10 -4 there is about 2 to 2.4 dB coding 
gain. However, there is an exception, the uncoded system responds better 
than the convolutional code when there is 100 % RFI. This result is due to 
typical limitations on the RFI energy. The magnitude of the RFI vector may 
not be sufficiently large to cause a bit error. 

Comparing the two TCM systems, the 16 State machine performs 
better than the pragmatic system (0.5 dB difference at 50% RFI). This 
difference in performance is due to the uncoded outboard bit used in the 
pragmatic system. It is, however, well worth noting that the pragmatic 
system is not optimal but, provides about 1.7 dB of coding gain compared 
to uncoded QPSK. 

4.7 BURST INTERFERENCE 

Perhaps the most popular source of man-made burst interference is 
a pulsed radar transmitter. Lightening is a natural phenomena also known 
to introduce random burst interferences. Therefore, communication 
transmissions over these channels introduces bursts as well as additive 
white Gaussian noise to message signals. Assuming that the occurrence of 
bursts are ideally discernible, they are treated as erasures in the received 
data. Parameters of the bursts that are of interest to the communicator are 
its duration and arrival interval. Two models of bursts interference were 
looked at. The first case modeled the burst as periodic pulses with regular 
arrivals. In this case the burst can be varied to control the length 'B' 
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BER 


(number of symbols erased) and the period 'Pb' in which it repeats. For the 



Composite Es/No dB 
Signal Power/(AWGN + RFI) 


Figure 4.6-1 Bit Error Rate vs SNR for QPSK. 
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Figure 4.6-2 Bit Error Rate vs SNR for 16-State Convolutional 

Encoder. 
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Figure 4.6-3. Bit Error Rate vs SNR fo Pragmatic System. 

second case the occurrence of a burst was random and the inter-arrival 
time had a Poisson distribution. The parameter 'B' as before controls the 
length while X controls the average arrival time. The general burst 
structure is shown in Figure 4.7-1. 



Burst length 
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Fi gure 4.7-1 Burst Parameters 
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Figure 4.7-2 RFI Burst Model 

BOSS modules were designed to generate both burst characteristics 
described above (see Figure 4.7-2). A pulse generator was used to produce 
the effects of recurring bursts. The pulse generator has output levels of 
zero and one. During the effected burst, the components of the signal 
vector are multiplied by zero or effectively erased. The rationale of 
designing the simulation in this fashion is as follows. True burst 
interference is of sufficient energy to render the information from the 
received signal vector essentially useless. This is the same effect as an 
erasure. Message symbols are erased when multiplied by the zero level of 
the pulse and this is done for the desired burst length. The burst 
interference sources were implemented with the TCM systems and BER 
were measured. During the occurrence of each burst 'B' symbols are 
erased. It is well known that convolutional codes responds poorly to these 
types of correlated errors 8 . Figures 4.7-3 and 4.7-4 show the performance 
curves both TCM systems with periodic bursts. Plots show, parametrically, 
the average number of incorrectly decoded output bits as a function of 
burst length at fixed signal-to-noise ratio. The degradation due to the burst 
is well illustrated by these plots. In low signal-to-noise ratios the bit errors 
are in excess of the number of erased symbols. However, as the SNR is 
increased the Viterbi decoder is able to correct erasures. This ability to 
correct erasures declines as the burst length increases. This error 
correction, even in the higher SNR, is an advantage compared to an 
uncoded QPSK system which offers no protection at all. 
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Figure 4.7-3. Bit Errors vs Burst Length for 16-State System. 
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Figure 4.7-4. Bit error vs. Burst Length, Pragmatic System. 
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For the uncoded QPSK system the probability of a bit error (with 
AWGN & periodic burst) is given by 


(1 - 5)Q 



+ 



where S is the duty cycle of the burst. For high values of Eb/N 0 the error 
probability is dominated by the burst duty cycle. If the arrival of the burst 
is random with a Poisson distribution, the error probability is obtained as 
follows. Let X be the average arrival rate of the burst. The probability of a 
single burst occurring in this time is given by 

p = 1- (1 - e-*) B , 


where B is the length of the burst. The 5 term in the previous equation is 
replaced by the expression ’p’, above, to give 


(1 - p) Q 



1 

+ 2 p * 


Once the QPSK system is subjected to burst interference its BER 
performance is limited by the nature of the burst. No amount of 
reshuffling or spreading (interleaving) of the burst will help. For a 10% 
duty cycle burst the bit error probability for QPSK approaches 0.05 
(independent of Eb/N 0 ). TCM systems, unlike QPSK, responds better when 
bursts are interleaved. 


4.8 PERIODIC CONVOLUTIONAL INTERLEAVING 

The periodic convolutional interleaver (PCI) used in this study was 
proposed by G. D. Forney. (See Figure 4.8-2). Convolutional interleavers 
offer some advantages over the more familiar block interleavers. Less 
memory storage, and shorter message delays are the major advantages of 
PCI over block interleaving. An interleaving example follows. Generally a 
communication system transmits sequential messages. If a burst is five 
symbols long, then each time a burst occurs five consecutive message 
symbols will be erased. As discussed previously, the Viterbi decoder will 
correct erasures provided that the interval between them is sufficiently 
long. The purpose of the interleaver is to scramble the transmitted 
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sequence so as to increase the distance between erased symbols at the 
receiver (see Figure 4.8-1). 


Burst Length 


Burst after interleaving/de-interleaving 


Distance Between 
Erased symbols 

Figure 4.8-1 Interleaver/De-interleaver Response to Bursts 

If the distance between erased symbols is sufficiently large the decoder is 
likely to perform much better than had the erasures remained clustered. 

Figure 4.8-2 shows the basic structure of the convolutional 
interleaver and its output. The de-interleaver at the receiver re-orders the 
received symbols to their original positions. Each line in the 

interleaver/de-interleaver has shift registers used to delay the input 
symbols. The first and last lines of the interleaver have no delay. 
Interleaving is accomplished by the action of the synchronized commutator 
and the various delays introduced by registers. 'M' is the smallest multiple 
of shift registers on any line. The interleaver depth *D’ is the maximum 
burst length that can be tolerated. N' = D * M is the minimum distance 
between consecutive symbols. The distance between consecutive erased 
symbols is however D * M - 1. This is illustrated in Figure 4.8-2 which 
shows the interleaved channel input. If a burst occurs, it affects the 
column of the matrix where each element in a column is separated by D * 
M - 1. 

In Figure 4.8-2 the interleaver shown has a depth of 5. Also M is 
equal to one; therefore, the distance between erased symbols is four. By 
increasing M, the number of shift registers, the distance between erasures 
is increased ( so are the delays). 
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A PCI, with D = 5, of the structure described, was implemented in 
BOSS and tested with the TCM systems. Transmissions were subjected to 
different levels of AWGN as well as a periodic burst with a 10% duty cycle 
(five symbols erased every 50 symbols). Bit error rates were measured for 
different values of N' and plotted versus Es/No- See Figures 4.9-1, 4.9-2, 
and 4.9-3. 
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Figure 4.8-2 Periodic Convolutional Interleaver 
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As before the abscissa Es/No is due only to additive white Gaussian noise 
and does not include any power due to bursts. In this paper interleavers 
are labeled as PCI(D,M) where D is depth and M is the smallest multiple of 
shift registers on a line. This was done to track the complexity of the 
interleaver. The BER measurements are plotted parametrically on this 
basis. 

4.9 BURST INTERFERENCE RESULTS 

For a (5,50) burst (10% duty cycle) five of every symbols are erased. 
The maximum distance that can be achieved between consecutive erased 
symbols is therefore ten (10). That is to say the (5,50) burst is transformed 
into a (1,10) burst. The higher order of interleavers (greater than 
(PCI(5,2)) provides distances greater than 10 within a single burst. 
However, because the bursts repeat periodically this distance is not 
maintained and can even revert to the non-interleaved case. It is to be 
expected that the best BER performance is achieved when the distance 
between erasures is 10 symbols. The BER curves obtained indeed support 
this conclusion (see Figure 4.9-1). Results were also obtained for random 
occurrence of bursts. Periodic interleavers become vulnerable to variations 
in the arrival rate and length of bursts and therefore perform worst when 
subjected to random bursts. Results also show that higher order 
interleavers provided better improvements for random bursts when 
compared to performances in periodic bursts. 

No amount of interleaving of the bursts recovered any losses for the 
Pragmatic TCM system. This result is due to the uncoded bit that is 
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Figure 4.9-1. 16-State TCM BER Performance with 10% duty 
Cycle Burst and Periodic Convolutional Interleaving 
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Figure 4.9-2. Bit Error Rate vs Es/NO for Burst with Poisson 

Distributed Arrival Rate 
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Figure 4.9-3 Pragmatic TCM BER 
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characteristic of the Pragmatic codes. The performance of the code is 
dominated by the outboard bits and interleaving does not provide any 
added protection (see Figure 4.9-3). 

4.10 CONCLUSIONS 

Trellis coded modulation (TCM) has been sighted as the modulation 
technique for the next generation of space communication systems 5 & 6. 
However, questions have surfaced concerning how well TCM systems will 
respond to radio frequency interference (RFI). 

This paper presents results of simulations geared to measure the 
performance of TCM systems in radio frequency (RF) co-channel 
interference and burst interference. The software package, Block Oriented 
Simulation Software (BOSS), was used to simulate (at rate 2/3) al6 State 
TCM system with Viterbi decoding and a Pragmatic TCM system. The 
performance of the system is measured by their bit error rates (BER) as a 
function of signal to effective interference ratio. 

Results show that when compared to the uncoded, quadrature phase 
shift keying (QPSK) system, TCM achieves a substantial coding gain when 
subjected to a mixture of AWGN and co-channel interference. This coding 
gain is in fact comparable to gains obtained with strictly Gaussian noise. It 
is a known fact that BER performance of TCM systems degrades 
significantly as a result of burst interference. However, periodic 
convolutional interleavers (PCI) can recover some of the losses. 
Interleaving is unsuccessful in recovering any losses when Pragmatic TCM 
systems are subjected to burst interferences. This result is due to the 
uncoded outboard bits. 

It can be stated that, under practical conditions, a properly designed 
Trellis Coded Modulation system will out-perform an uncoded QPSK system 
in environments subjected to both Gaussian and co-channel interference, 
and burst interference. 
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SECTION 5. CONCATENATED TRELLIS-CODED MODULATION 
5.1 INTRODUCTION 

Concatenated systems combining the error correcting power of trellis 
code and Reed-Solomon (RS) codes have been used to achieve error rates 
between 10' 10 and 10'20 on satelite communications links [1]. Typically 
such systems use binary or quaternary signalling, with an industry 
standard binary convolutional codes. The need for spectral efficiency has 
motivated consideration of of non-binary signalling with Ungerboeck [2,3] 
codes or the Pragmatic standard recently proposed by Viterbi [4]. In this 
paper we investigate the use of 8-PSK trellis-coded modulation (TCM) in a 
concatenated system. The outer code will be the (255,233) RS code, as this 
is the standard for current systems. This investigation concentrates on the 
effect of substituting TCM for the currently used binary inner code. 



Figure 5.1-1. Concatenated Code System 


A concatenated system is shown in figure 5.1-1. User data is first 
encoded by the RS Block encoder. The RS code symbols are then 
interleaved and convolutionally encoded. At the receiving end of the 
channel, modulated signals are decoded using the Viterbi algorithm, 
preferably with soft decisions. The output from the Viterbi decoder is 
then deinterleaved and RS decoded, recovering the user data. The RS code 
block consists of 255 symbols, a symbol being a block of 8-bits. An error 
in one or more of the eight bits of a given symbol results in a symbol error. 
The RS block code can correct up to 11 symbol errors. 

The purpose of the interleaver is to aleviate the effect of burst noise 
during transmission. Given the structure of the RS code, we can see what 
the proper operation of the interleaver should be. Since a single bit error 
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in an RS code symbol counts as a symbol error, additional bit errors in the 
same symbol cause no further damage. An interleaver-deinterleaver 
operation which redistributes bit errors among multiple symbols in the 
same block will degrade, not improve, the performance of the system. 
Therefore, the interleaver should be designed to reduce the correlation 
between symbol errors, and to redistribute clustered symbol errors among 
several RS blocks. This fact will be important in the forgoing analysis. 

Needless to say, it is not practical to measure bit error rates between 
10-1° and 10' 20 by computer simulation. When the concatenated system 
is operating in this range, the Viterbi decoder will be operating at a bit 
error rate between 10' 3 and lO" 5 . Therefore, it is feasible to characterize 
the Viterbi decoder output using simulations, and then algebraicly 
determine the resulting probability of error for the RS code. Such an 
approach has been used to predict the performance of concatenated 
systems with binary signalling [5], here we extend the technique to 
concatenated systems using TCM. 

5.2 CHARACTERIZATION OF THE VITERBI DECODER OUTPUT 

The approach used in [5] suggests that the output of the Viterbi 
decoder be modelled as a two-state first-order Markov process. The 
process, illustrated by the state diagram of figure 5.2-1, consists of a burst 
state, in which the probability of error is very high, and a clear state, in 
which the probability of error is very low. The process is characterized by 
the probabilities of transition from either state to the other. The 
transitional probabilities reflect the intercorrelation between errors in the 
Viterbi decoder output. 


F*bc 



Figure 5.2-1. State Diagram for Markov Process Model 

The Markov process is completely specefied by four parameters: 

P|b = probability of error when in burst state (5.2-1) 

Pic = probability of error when in clear state (5.2-1) 
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Pbc = probability of transition to clear state, 
when in burst state 

P C b = probability of transition to burst state, 
when in clear state. 

(5.2-3) 

(5.2-4) 

From these four, additional useful parameters may be 

derived: 

Pgb = I-P15, the probability of non-error, 
when in burst state 

(5.2-5) 

Pq c = 1-Pi c , the probability of non-error, 
when in clear state 

(5.2-6) 

Pbb = 1-Pbc> ^e probability of remaining 
in burst state 

(5.2-7) 

P cc = 1-P c b, the probability of remaining 
in clear state 

(5.2-8) 

Pc b 

Pb = D D > the probability that the process 

*cb + *bc 

(5.2-9) 

will be in the burst state at any given time 


Pbc 

P c = D d » the probability that the process 

Pcb + *bc 

(5.2-10) 

will be in the clear state at any given time. 


To determine these parameters, which depend on channel signal-to- 
noise ratio, it is necessary to observe statistics of bit error patterns in the 
Viterbi decoder output. We define the following observable 

characteristics: 

Pj = the probability of a single bit error 

(5.2-11) 

P 1 j = the probability of two consecutive bit 
errors 

(5.2-12) 

Pill = the probability of three consecutive bit 
errors 

(5.2-13) 

PlOi = the probability of a bit-error, a 

(5.2-14) 
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non-bit-error, and a bit-error, in 
succession. 

The observable characteristics are algebraicly related to the process 
parameters as follows: 

Pl= Pb p lb + p c p lc (5.2-16) 

Pi 1 = p bPlb p bb p lb + Pb p lb p bc p lc (5.2-17) 

+ P c PlcPcbPlb + PcPlcPccPlc 

Pill = Pb p lbPbbPlb p bb p lb + p bPlb p bcPlc p cbPlb (5.2-18) 

+ ... + P c PlcPccPlcPcc p lc 

p 101 = p b p 1 b p bb p ObPbbP 1 b + PbPlb p bc p OcPcbPl b (5.2-19) 

+ ... + PcPlcPccPOcPccPlc 

The expressions for Pi i and p i 01 are similar to equations (5.2-16) and 
(5.2-17) except that they consist of eight terms of eight factors each. 

If the two-state process is a valid model of the Viterbi decoder 
output, then the expressions for Pi, p ll» P 1 1 1 and p 101 (5.2-16 through 

5 2-19) are four equations from which the four unknowns Plb, Pic, Pbc, 
and P c b may be determined. The system is nonlinear, and there will in 
fact be two solutions, one of which exactly reverses the roles of the burst 
state and the clear state. Also, if p bc = p cc and ^cb = p bb> or p lb - lc> 
the model will generate a process in which bit errors are independent, and 
the solution will not be unique. 

Given a set of numbers Pj, p i l , p l 1 1 and P 1 0 1 not actually 
generated by the two-state process, a real solution may or may not exist. 
The algebraic solution of the system (5.2-16) through (5.2-19) is an 
exceedingly difficult exercise, and a reasonable approach to finding the 
model parameters is to use a numerical search to find the process 
parameters which match the observed characteristics. This approach was 
used to determine the parameters which charaterize the output of a TCM 

decoder. 

5.3. CALCULATION OF THE PROBABILITY OF SYMBOL ERROR 

Given the parameters of the Markov process, it is possible to 
calculate P(B), the probability that the process will be in the burst state B 
times during a symbol conisting of K > 3 bits. A symbol error is defined to 
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be the event that one or more of the bits of the symbol are in error 
Therefore, the probability of symbol error is given by 


p s = X( 1 - p 0b Bp 0c K - B ) 

B=0 


(5.3-1) 


P(B) is found by considering the K-bit symbol to be divided into 
segments of consecutive bits in burst state or in clear state. If we let Rb he 
the number of segments in burst state, and R c be the number of segments 
in clear state, then it is clear that 0 < Rb < B and R c = Rb - Rb or Rb + 1. 


Now, 


P(B,Rb) = P(B,R c =Rb-l) + P(B,Rc=Rb) + P(B,Rc- R b+l> (5.3-2) 

(5.3-3) 


P(B,R c =Rb-D 

P(B,R b =R c ) 

( B-l YK-B-1V R b -1 D Rc d 

l^Rb^ARc-i / bPcb Pbc Pbb 


R b-1 D . R c D , , B-Rbp K-B-R C 


cc 


(5.3-4) 


Rb-1 D . Rc d , . B-Rbp K-B-Rc 


CC 


( B-l YK-B-n Rb n Rc-1 d B R bp R -B R c 
R b -l) R c -l fcPcb ^ Pbb Pec 


P(B,R b +D v 

(K-lYK-n n R bo R c -1. B-Rbp K-B-Rc 
^Rb-ljRc-l/^h Pbc Pbb CC 

The reasoning behind this is as follows. The combinatorial 
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Rb-^j is the number of wa y s that B p eriods in the burst state can be 


(K-B-n 

divided among Rb segments in burst state, and likewise ^R c -1 


is the 


number of ways that K-B periods in clear state can be divided among R c 
segments in clear state. Every segment except the first segment begins 
with a transition, hence the factor P c b or Pbc is used - The first se 8 ment 
begins with Pb or P c . The remaining K periods must be non-transitions, 
hence P b b or P cc . Equations (5.3-1) through (5.3-5) are valid for 1 < B < K, 
Rb > 1 and R c > 1. If Rb=0 then B=0 and if R c =0, B=K. 


P(B=0) = PcPcc 1 ^' 1 


(5.3-6) 


P(B=K) = PbPbb K ’ 1 


(5.3-7) 


ForO<B<K: 


K-l 

P(B) = ^ P(B, R b) 

R b =l 


(5.3-8) 


Where P(B,Rb) is as specified by equations (5.3-1) through (5.3-5). 
Once P(B) is calculated using equation (5.3-8), the probability of symbol 
error is calculated using equation (5.3-1). 

5.4. CALCULATION OF THE REED-SOLOMON BLOCK ERROR 
PROBABILITY 

The (255,233) Reed-Solomon code consists of a block of 255 symbols 
of 8-bits each’. Up to 11 symbol errors are correctable therefore, the 
probability of Block error is given by 

Pe-l-£ (T^O-Ps) 255 -* (5 ' 4 -' ) 

For completely independent bit errors: 

P s = 1 - Pi 8 (5.4-2) 
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For bit errors described by the model of section II, P s is calculated as 

described in section III, using K=8. 

If the RS codebits are transmitted using a rate 2/3 8-PSK Ungerboec 
code each state transition of the Markov process must correspond to a 
period of operation of the trellis code, which corresonds to two decoded 
bits. Therefore, we introduce the term Viterbi decoded symbol to refer to 
a pair of bits which are decoded in one stage of the Viterbi decoder trellis, 
and a RS code symbol then consists of four Viterbi decoded symbols. n 
determining the parameters from the Markov process model, a Viterbi 
decoded symbol is" conidered to be in error if either of its two constituent 
bits are in error, and this probability becomes the statistic Pi- Likewise, 
Pjl denotes the probability of two consecutive Viterbi decoded symbol 
errors, and Pill denotes the probability of three such errors. The 

parameters of the Markov process are then evaluated using the method of 
section II, and P s is evaluated as described in section III, using K-4. Ibe 
validity of equation (5.4-1) depends on the independence of RS code 
symbol errors, which is effected by the interleaver designed to interleave 

blocks of four 8-PSK symbols. 

5.5. RESULTS 

Simulations were written in the C programming language to 
determine the probability of all patterns of errors in four consecuttve 
Viterbi decoded symbols. The results for the 16-state Ungerboeck code 
are complete. The results for the pragmatic code and the 64-state 
Ungerboeck code are in progress. These results were used to calculate the 
probability of a block error in the Reed-Solomon code, wh.ch is plotted in 

figure 5.5-1. 
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Figure 5.5-1. RS Block Error Probability for Concatenayed System 
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